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The term of avalanche may originate from the 
phenomena occurred in nature. It is referred to 
as a sequential events which may cause devastat- 
ing catastrophe. The phenomena of avalanches 
are ubiquitous in nature. The canonical exam- 
ple of avalanche in nature is the mountain slide, 
during which great mass of snow and ice at a 
high altitude slide down a mountain side, often 
carrying with it thousands of tons of rock, and 
sometimes destroy forests, houses, etc in its path 

[I] . Since avalanches occur everywhere, from the 
sandpile or ricepile, to the Himalayan sandpiles; 
from the river network, to the earthquake, star- 
quakes and even solar flares; from the biology, to 
the economy, [2] , etc, it is hence proposed [2] that 
avalanches may be the underlying mechanism of 
the formation of various geographical structures 
and complex organisms, e.g., brains, etc. Further- 
more, avalanches may be the origin of fractals in 
the world. From this point of view, avalanches 
can be viewed as the immediate results of complex 
systems, and hence can be used as the theoretical 
justification for catastrophism. This is because if 
the real world is complex then the catastrophes 
are inevitable and unavoidable in biology, history 
and economics. It is now even proposed by Meng 
et al [3] that the formation of colorless gluon clus- 
ters may be attributed to avalanches intrigued by 
the emission or absorption of gluons. 

Plenty of patterns provided by nature exhibit 
coherent macroscopic structures developed at var- 
ious scales and do not exhibit elementary inter- 
connections. They immediately suggest seeking 
a compact description of the spatio-temporal dy- 
namics based on the relationship among macro- 
scopic elements rather than lingering on their in- 
ner structure [4]. That is, one needs to condense 
information when dealing with complex systems. 
Maybe only this way is efficient and turns out suc- 
cessful. 

As known, avalanche is one kind of macroscopic 
phenomenon driven by local interactions. The size 
of avalanche, spatial and temporal as well, may be 
sensitive to the initial configuration, or more gen- 
erally, the detailed dynamics of the system. How- 
ever, the distribution of avalanches, Gutenberg- 
Ritcher law [5], or equivalently, power-law, does 
not depend on such kind of details due to the 
universality of complexity. Hence, in this sense 
avalanche study may be an appropriate tools in 
studying various complex phenomena. On the 
other hand, observation of a great variety of pat- 
terns, such as self-similar, fractal behavior in na- 
ture [6,7,8,9], J noise in quasar [10], river flow 

[II] and brain activity [12], and many natural 



and social phenomena, including earthquakes, eco- 
nomic activity and biological evolution suggests 
that these phenomena are signatures of spatiotem- 
poral complexity and can be related via scaling re- 
lations to the fractal properties of the avalanches 
[13]. This suggests the occurrence of these gen- 
eral, empirical phenomena may be attributed to 
the same underlying avalanche dynamics. Thus, 
one can see that study of avalanche is crucial in 
investigating the critical features of complex sys- 
tems. It can be even inferred that avalanche dy- 
namics docs provide much useful information for 
us to understand the general features of the ubiq- 
uitous complexity around us. That is probably 
also why this paper focuses on such kind of topic. 

Despite the fact that avalanche may provide in- 
sight into complexity, the definition of which can 
be vastly different for various systems, and the 
same kinds of systems, even the same system. Let 
us recall some definitions of avalanches given be- 
fore. In sandpile model [2], an avalanche is in- 
trigued by adding a grain or several grains of sand 
into the system at some time and causing the top- 
ple of some sites, which may later on cause some 
other sites to topple. The avalanche is considered 
over when the height of all sites are less than the 
critical value, say, 4. In Bak-Sneppen model [14], 
several kinds of avalanche [13] are presented. For 
instance, /o-avalanche, G(s)-avalanche, forward 
avalanche, backward avalanche, etc. Despite the 
fact that these kinds of definitions of avalanches 
may show the various hierarchical structures they 
manifest the same underlying fractal feature of the 
ecosystem, i.e., self-organized criticality. Relating 
all these kinds of avalanches one can provide a gen- 
eral definition of the avalanche for Bak-Sneppen 
model: An avalanche corresponds to sequential 
mutations below certain threshold. One can see 
that this kind of definition can ensure the muta- 
tion events within a single avalanche are casually 
and spatially connected. In addition, with this 
definition there exists a hierarchy of avalanches, 
each defined by their respective threshold. It is 
the hierarchical structure of the avalanche that ex- 
hibits the fractal geometry of the system and that 
implies complexity. 

It can be inferred from the definition of 
avalanche that there always exists a triggering 
event which initiates the avalanche and whose ef- 
fect, that is, causing an avalanche to spread within 
the system later and later on, will disappear at 
the end of the avalanche. And, the observation 
of avalanche through the triggering event, up to 
now, is based on the individual level, despite that 
the avalanche is a macroscopic and global phe- 
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nomcnon of the system studied, in the laboratory, 
and in nature as welL Take sandpile model, the 
triggering event is adding a grain or several grains 
of sand to some sites and causing them to top- 
ple, thus initiating an avalanche. Consider an- 
other model, Bak-sneppen model, in which the 
corresponding triggering event of an avalanche is 
mutation of the extremal species causing the fit- 
ness [14] of the extremal site at the next time step 
less than a certain threshold. One can see that in 
the above two models triggering events are directly 
concerned with the feature of individuals, e.g., the 
height of the site in the former model, or, the fit- 
ness of the extremal site in the latter one. It can be 
readily learned that the triggering events, whether 
those in the laboratory or those in nature, are not 
directly related to the global feature of the sys- 
tems although avalanche can span across the whole 
systems. Generally speaking, the observation of 
avalanches is done through some feature of indi- 
viduals, instead of that of the system as a whole. 
However, general feature of the complex system 
may provide insight into knowing the tendency of 
the evolution of the system. Specifically, global 
feature of a complex system may help one to un- 
derstand the critical behavior of the system. That 
is, it is feasible that some characteristic quantities 
representing the corresponding global features can 
be employed in describing the critical behavior of 
the system. Furthermore, these quantities ought 
to be related to avalanche dynamics, and hence 
can be used to describe complexity emerged in a 
variety of complex systems. Apparently, our aim is 
to search for or define such kind of quantities and 
then to expect to observe new kind of avalanche 
based on these quantities. Indeed, we obtain a new 
quantity which can be used to define a new hier- 
archy of avalanches in Bak-Sneppen model. We 
suggest that this quantity may be used as a crite- 
rion in determining the emergence of criticality. It 
will be shown later that this new kind of avalanche 
still exhibits spatio-temporal complexity in a dif- 
ferent context. 

Consider Bak-Sneppen model [14], which is a 
very simple evolution model of biological ecosys- 
tem. Despite the simplicity of the model itself, 
it can exhibit the skeleton of species evolution, 
punctuated equilibrium behavior. Detailed infor- 
mation about this original model of evolution can 
be available in Ref. [14]. In Bak-Sneppen model, 
each species is represented by a single fitness. The 
fitness may represent population of a whole species 
or living capability of the species [15]. Hence, one 
can see that fitness is a vital quantity and is the 
only one describing the model. No other addi- 



tional quantities are considered in this oversimpli- 
fied model. Thus, the fitness is the most impor- 
tant feature of species and that of the model. So, 
when considering the global feature of the species 
ecosystem, one has to relate this general feature 
to the feature of individuals. That is, the gen- 
eral feature of the ecosystem should be associated 
with the fitness of the species. As previously men- 
tioned, a corresponding quantity should be found 
to describe this general feature. Before presenting 
such kind of quantity let us briefly review Bak- 
Sneppen model so that the readers who are not so 
familiar with this model can have a rough idea of 
what it is about. 

Bak-Sneppen model is perhaps the simplest 
model of self-organized criticality. In this "toy" 
model, random numbers, fi, chosen from a fiat 
distribution, p{f), are assigned independently to 
each species located on a d-dimensional lattice of 
linear size L. At each time step, the extremal site, 
i.e., the species with the smallest random number, 
together with its 2d nearest neighboring sites, is 
chosen for updating by assigning 2d + 1 new ran- 
dom numbers also chosen from the same uniform 
distribution to them. This updating process 
continues indefinitely. After a long transient pro- 
cess the system reaches a statistically stationary 
state where the density of random numbers in the 
system vanishes for f < fc and is uniform above 
fc, the self-organized threshold. 

Having briefiy introduced the model, next, we 
will introduce a new quantity. Please note that 
the model we used is still Bak-Sneppen model. We 
observe the evolution of the model without adding 
anything to the model. We simply introduce the 
new quantity based on the fitness of the species. 
Define the average fitness, denoted by /, as below, 

i=l 

, where fi is the fitness of the ith species. Here, 
we refer to / as the average fitness of the whole 
system and as a global quantity. / may represent 
the average population or average living capability 
of the whole ecosystem. Large /, i.e.. high aver- 
age fitness, may imply the total population of the 
system is immense or its average living capability 
is great, and vice versa. Initial value of /, denoted 
by /(O) , can be easily calculated. As known, at 
the beginning of the evolution /j's are uniformly 
distributed between (0,1). So, for an infinite sys- 
tem, /(O) equals to 0.5. However, for a finite-size 
system /(G) will fiuctuate slightly due to the fi- 
nite size of the system, which is not so important 
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in the latter evolution. We will simply consider 
the average value, 0.5. It should be pointed out 
that /(O) does not reflect the correlation among 
species. As the evolution goes on the correlation 
among the species within a system will become 
more and more distinctive. Denote f(s) the av- 
erage fitness of the system at time step s in _the 
evolution. Hence, in the s limit, i.e., s » L"^, f{s) 
may partly reflect information about correlation. 
As a global quantity, /(s) should include informa- 
tion concerning the interaction between species. 
Hence, it is natural to expect that / may be a 
good quantity in describing the feature of the sys- 
tem as a whole. 

Before introducing the new hierarchy of 
avalanches it is necessary and worthwhile to know 
some feature of the new quantity, /(s). Firstly, 
let us present some theoretical analysis. Recall 
the definition of / one can see that A/(s) = 
f{s + 1) — f{s) approaches zero in the L oo. 
An observer can not even perceived the change of 
f{s) during the short time period since it is van- 
ishingly small . However, changes at very time 
step are accumulated to form a relatively distinc- 
tive change after a long time, which is perceiv- 
able for the observers. This long time period is 
required to be much greater than the system size, 
i.e., s :$> L"^. In other words, f{s + so) — f{so) may 
only be "noticed" when s L'^. Thus, one can 
not expect f{s) will have great variation from the 
current time step to the next nearest time step, 
which is vastly different from the variation of the 
fitness of extremal site. The latter can vary from 
one value, say, 0, to the next value, 1 between two 
successive time steps. It should also be expected 
that there exists an increasing tendcy of / ver- 
sus time s. This is because at each time step the 
least fitness is eliminated from the system so the 
general fitness of the whole system will tend to in- 
crease. And due to the slow fluctuation of the / 
the increasing behaves like a stepwise, i.e.. Devil's 
stepwise [2] . One then may expect to observe such 
behavior, i.e., punctuated equilibrium [14], of / in 
the evolution of Bak-Sneppen model. 

In order to show the feature of / versus time s 
we performed simulations of Bak-Sneppen model. 
At each time step, in addition to the updating of 
the extremal sites, we also track the signals J{s). 
FiG. 1 shows the evolution of /(s) versus time 
during a time period for a one-dimensional Bak- 
Sneppen model of size L = 200. This plot shows 
that / varies slightly between two successive time 
steps but tends to increase in the long evolution 
process. Simulation of a two-dimensional model 
of size L = 20 exhibits the similar behavior of the 



evolution of f{s). 

Before searching for the punctuated equilibrium 
behavior let us flrst introduce another quantity, 
F(s), the gap of the average fitness. The defini- 
tion of F(s) is given as follows: Initial value of 
F(s) is equal to /(O). After s updates, a large 
F{s) > F{0) opens up. The current gap F(s) is 
the maximum of all F{s'), for all < s' < s. FiG. 
2 shows the F(s) as a step-wising increasing func- 
tion of s during the transient for a one-dimensional 
Bak-Sneppen model of size L = 100. Actually, the 
gap is an envelope function that tracks the increas- 
ing peaks in /(s). Indeed, punctuated equilibrium 
behavior appears in terms of this new quantity, 
lis). 

By definition [14], the separate instances when 
the gap F(s) jumps to its next higher value are 
separated by avalanches. Avalanches correspond 
to plateaus in F(s) during which /(s) < F{s), 
which ensures the mutation events within a single 
avalanche are spatially and casually connected. A 
new avalanche is initiated each time the gap jumps 
and ends up when the gap jumps again. As the 
gap increases, the probability for the average fit- 
ness, /, to fall below the gap increases also, and 
larger and larger avalanches typically occur. 

We can also obtain an exact gap equation of 
F(s), similar to the one found for Bak-sneppen 
model in Ref. [16]. Suppose in the system the cur- 
rent gap is F(s). If F(s) is to be increased by AF, 
i.e., from F(s) to F{s) + AF, the average number 
of avalanches needed is A^av = AFL'^/{1 — F{s)). 
We can guarantee A^av ^ 1 by selecting AF ^ 
L~'^. In the large L limit, A^av can be arbitrar- 
ily large. Hence, in this limit, the average num- 
ber of time steps required to increase the gap 
from F(s) to F{s) + AF is given by the inter- 

Val As = (5)F(s)^av - (s}F(s)AFiV(l - F{s)), 

where (S')p(s) is the average avalanche size of the 
plateaus in the gap function. From the law of large 
numbers the fluctuation of this interval around its 
average value vanishes. In the AF limit. 
As 0. Taking the continuum limit we can ob- 
tain the differential equation for F(s), 

dF(s) ^ 1 - F(s) 
ds i'^(5>F(,)- ^ ' 

Note this equation is exact. 

All SOC models, e.g., the BTW sandpilc model 
[17], the earthquake models [18], or Bak-Sneppen 
model [14], exhibit self-organized criticality in 
terms of a power-law distribution of avalanche. It 
is natural to expect that we can observe SOC in 
terms of the hierarchical structure of /, which it- 
self manifests complexity. Using this new quantity 
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to define the avalanche is simply another way of 
observing the same phenomenon which can be ob- 
served in other ways. As known, the emergence 
of complexity is independent of the tools used to 
observe them provided that these tools are effi- 
cient and strong enough. Similar to the ones used 
in Refs. [13,19], we present the definition of Jo- 
avalanche, where /o is only a parameter between 
0.5 and 1 to define the avalanche. Suppose at 
time step si, /(si) is larger than /q. If , at time 
step si-l-1, /(si -|- 1) is less than /o, this initiates 
a creation-annihilation branching process. The 
avalanche still continues at time step s', if all the 
/(s) are less than /o for 1 < s < s' — 1. And the 
avalanche stops, say, at time step Si -|- S, when 
/(si + S) > f{si). In terms of this definition, 
the size of the avalanche is the number of time 
steps between subsequent punctuation of the bar- 
rier /o by the signal /(s). In the above example, 
the size of the avalanche is S. It can be clearly 
seen from FiG.l that this definition guarantees 
the hierarchical structure of avalanches, larger 
avalanches consists of smaller avalanches. As /q 
is lowered, bigger avalanches are subdivided into 
smaller ones. Hence, the statistics of /o-avalanche 
will inevitably have a cutoff if /o is not chosen to 
be the value of /(s) at critical state, denoted by fc- 
We can also define the /c-avalanche. Nevertheless, 
/o-avalanche in the stationary state has the same 
scaling behavior as /c-avalanche provided that /o 
is close to fc- We measure /o-avalanchc distribu- 
tion for one- and two-dimensional Bak-Sneppen 
models. The simulation results are given in FiG. 
3. The exponent r, defined by P{S) ^ s~'^, is 
1.80 for ID model and 1.725 for 2D model. An- 
other exponent, D, avalanche dimension [13], de- 
fined by ricov ~ S'-'/'^, where ricov is the number of 
sites covered by an avalanche, and d is the space 
dimension, is measured. We find £'=2.45 for ID 
model and 1.55 for 2D model. 

Up to now, a question is_ still unsolved. It is 
about the critical value of /, fc- This may be a 
hard bone if the system size is finite, but when the 
consider the L limit, everything will be smooth 
and can be easily accomplished. Recall the evo- 
lution of Bak-Sneppen model, or the detailed re- 
search of this model [13], the densities of sites with 
random numbers is uniform above G [13] and van- 
ishes below G in the L oo, where G is the gap 
of extremal site and detailed information of it can 
be found, for instance, in Ref. [13]. Hence, one 
can obtain, 



tion of G found in Ref. [19], one can immediately 
obtain Eq. (2). Please note that Eq. (2) is also 
valid for finite-size systems. From Eq. (3) one can 
immediately obtain, 

l + /c 



lim /c = lim 



(4) 



lim /(s) = lim 



l + G{s) 



(3) 



Interestingly, inserting Eq. (3) into the gap equa- 



Hence, fc can be easily determined from Eq. (4). 

Using the results of fc provided by Refs. [13,20], 
one can obtain /c, 0.83351 for ID model and 
0.66443 for 2D model. However, Eqs. (3) and 
(4) are not valid for a finite-size system, since one 
can not ensure the distribution of random numbers 
during a finite-size system is really uniform. Due 
to the fiuctuation of / (s) it is extremely difficult 
for one to determine exactly the critical value of / 
for a finite-size system. One may estimate fc for 
a finite-size system using the simulation. We find 
that this value weakly depends on the system size. 
When the system size is very large fc approaches 
the corresponding value for infinite systems. Actu- 
ally, the value of fc itself is not so important. FiG. 
4 shows the fluctuation of / for a one-dimensional 
model of size L = 200 near its critical state. We 
note, in this figure, / fluctuates slightly around 
some average value and does not tend to increase 
any more for a long time. We may say that the sys- 
tem approaches its stationary state. In this sense, 
we suggest that / may be a good quantity in de- 
termining the emergence of criticality. That is, the 
great fluctuation of /min will not affect us to de- 
termine when we approach the critical state. We 
need only to know the feature of /. This is more 
reasonable and easily accepted since / is a global 
quantity and condenses information of the system 
and its components. 

Why we call the /o-avalanche a new hierarchy 
of avalanches? Firstly, this kind of avalanche is 
deflned on the global level, in terms of the new 
global quantity, /. The background of this defi- 
nition is different from any one used before. This 
new kind of avalanche reflects the fractal geome- 
try in terms of the global feature. Secondly, one 
can notice that the exponents r obtained in our 
simulation are different from the ones found in 
Ref. [13]. From this point of view, one can judge 
that this kind of avalanche is totally different from 
any one observed before. Hence, its a new kind of 
avalanche. 

Self-organized criticality is suggested by Bak et 
al to be the " flngerprint" of a large variety of com- 
plex system (they call system with variability as 
complex) and is represented by a scale-free line 
on a log-log plot. In order to know the critical- 
ity of a system one needs to know when the sys- 
tem reaches the stable stationary state where the 
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phase transition occurs. It is extremely difficult 
and almost impossible for one to know when a sys- 
tem in nature approaches, not even reaches , its 
critical state. One can just study the ubiquitous 
fractal geometrical structure carved by avalanches 
through thousands of millions of years. However, 
in laboratory experiments and computer simula- 
tions, one needs a criterion to judge when station- 
ary state approaches, even, reaches, since statis- 
tics of avalanches may only be done under critical 
state of the system. Given Bak-Sneppen model, 
when the extremal signal, /mim approaches to the 
self-organized threshold, /c, the ecosystem reaches 
its stationary state. However, itself fluctuates 
greatly time to time, which brings great difficulty 
in determining the appearance of criticality. Thus, 
we provide a new quantity, /, for a candidate in 
judging the emergence of criticality. As shown, / 
is relatively stable in a short time period. Hence, 
when / does not tend to increase any more, we 
may say that the system approaches its station- 
ary state. And, we can observe criticality in a 
rather long time period. Surely, the emergence of 
criticality is rather complex, other physical mech- 
anism is needed, this is what we will consider in 
the future work. 

In conclusion, a new hierarchy of avalanches is 
observed in Bak-Sneppen model. A new quantity, 
/, is presented and is suggested by us to be a pos- 
sible candidate in determining the emergence of 
criticality. An exact gap equation and simulation 
results are also given. 
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FiG. 1: The variation of / versus time dur- 
ing a time period for a (a) one-dimensional Bak- 
Sneppen model of size L — 200 and (b) two- 
dimensional Bak-Sneppen model of size L — 20. 
The plots show the hierarchical structure of /. 

FiG. 2: Punctuated equilibrium of / for a (a) 
one-dimensional Bak-Sneppen model of size L = 
200 and (b) two-dimensional Bak-Sneppen model 
of size L — 20. We track the increasing signal of 
i.e, F(s). 

FiG. 3: Distribution of /o-avalanche for a (a) 
one-dimensional Bak-Sncppcn model of size L = 
200 and (b) two-dimensional Bak-Sneppen model 
of size L = 20. /o for (a) is chosen to be 0.821, and 
for (b), 0.648. The slopes are -1.800 and -1.725 for 
the two plots respectively. 

FiG. 4: The fluctuation of / around the critical 
state of a one-dimensional Bak-Sneppen model of 
size L = 200. 
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